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Abstract 

The modifications induced in the calculation of the cross section of the diffrac- 
tive process 77 — > J/^J/ip when the gluon propagator is changed are ana- 
lyzed. Instead of the usual perturbative gluon propagator, alternative forms 
obtained using non-perturbative methods like Dyson-Schwinger equations are 
used to consider in a more consistent way the contributions of the infrared 
region. The result shows a reduction in the differential cross-section for low 
momentum transfer once compared with the perturbative result, to be con- 
firmed with future experimental results from TESLA. 
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The description of the diffractive processes in high energy physics is still an open question. 
In the high energy limit, those processes are usually described by a Pomeron exchange [|TJ, 
whose nature remains one of the most intriguing phenomena in Quantum Chromodynamics 
(QCD). Its properties were first described by the Regge theory and used to explain the 
features of total hadronic cross-sections and differential cross-sections at small transverse 
momenta (the soft Pomeron), whose properties relied on the non-perturbative sector of the 
QCD. In high energy QCD, the Pomeron corresponds to the sum of gluonic ladder diagrams, 
with reggeized gluons in the vertical lines, described by the Balitskii-Fadin-Kuraev-Lipatov 
(BFKL) equation (the hard Pomeron) ||. Although the BFKL equation is infrared finite, 
there are difficulties to control contributions from the infrared region, reflecting in the fact 
that the BFKL approach does not give the same results as the Regge theory for the soft 
Pomeron. Since perturbative QCD can only be applied when all momenta are sufficiently 
large [J], in a region of low momentum transfer we can expect an overlap of the effects from 
the soft and hard sectors. 

In an experimental procedure to test the hard Pomeron, the observables should have 
some particular properties: the virtuality of the gluons in the ladders must be sufficiently 
large to enable the use of a perturbative expansion and the hard scale is provided by the 
coupling of the ladder to the incoming particles or by the momentum of the gluons in the 
ladder. 

Due to the requirements above, measurements in hadron-hadron and lepton-hadron col- 
liders are limited. Kwiecihski and Motyka j| argued that these difficulties can be avoided in 
the case of the production of massive vector mesons from two photon reactions at TESLA || . 
The exclusive process 77 — > J/^ J/^ (see fig. [0) can test the exchange of the Pomeron 
in arbitrary momentum transfer since the hard scale is given by the (large) mass of the 
quark charm, for both sides of the diagram, instead of the gluon momenta in the ladder. 
However, in the infrared region (low momentum transfer) where the soft Pomeron is more 
important, we expect non-perturbative effects will play a crucial role in the description of 
these processes. To estimate these soft effects, in this letter we propose as a first step the 
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replacement of the usual perturbative gluon propagator by distinct ones that include non- 
perturbative effects, considering the soft Pomeron as the exchange of two non-perturbative 
gluons. This model for the Pomeron was used successfully to describe another soft processes 
as the elastic proton-proton scattering using the gluon propagator from Dyson-Schwinger 
equations (DSE) |j] or from lattice field theory (LFT) |7J and the elastic electro-production 
of vector mesons |§. The modified propagators can be obtained through different ways: 
DSE H and numerical simulations in lattice field theory ||10|| , for example. 

The paper is organized as follows: first, we present the model to calculate the 77 — > 
J/^J/ty process using the approach of the ref. ||. Next, the main features of the employed 
non-perturbative gluon propagators are reviewed and then, we perform the calculation of 
the differential cross section, da/dt in the two gluon approximation with a modified gluon 
propagator, and finally present the conclusions and discuss some possible future develop- 
ments. 

The differential cross section for the 77 — ► J/^J/^ obtained in the BFKL approach, is 
given by @| 

where 

c> M *\- f d " k $o(fc 2 ,g 2 )$(x,k,Q) 

j J [(k+Q/2)2 + So ][(k-Q/2)2 + So ] [Z) 

and k ± Q/2 is the transverse momentum of the exchanged gluons; s = W 2 is the total 

center of mass energy of the 77 system; t = —Q 2 is the square of the transverse part of 

the momentum transfer (see figure [I]); x — m 2 j^/W 2 , with rrij/q the J/^f mass and s is 

a parameter related to the employed gluon propagator (we will discuss about this point 

below) . 

The impact factor for the transition 7 J/^ is $o(& 2 ; Q 2 ) which is induced by two gluons 
(see figure ||]) and in a non-relativistic approximation Jl2| is 



%{k 2 ,Q 2 ) = -^a^asifj 2 ) 



q 2 m 2 ji^/A + k 2 



(3) 



In the above formula 



C = gc-vrmj/vt/j/vK, (4) 

where qc = 2/3 is the charm quark charge; fj/y is a parameter that characterizes the 
light-cone wave function of the charmonium and considered here with value 0.38 GeV [Q; 
a em is the electromagnetic coupling constant; a>s(^ 2 ) = 127r/(25 ln(/i 2 /A 2 )) is the strong 
running coupling constant in the one loop approximation (the QCD scale is A = 0.23 GeV), 
f = {m 2 j^ + Q 2 )/4 and /x 2 = k 2 + Q 2 /4 + (mj/^/2) 2 . 

The function $(x,k, Q 2 ) obeys the BFKL equation in the leading ln(l/x) approxima- 
tion 01 



3a(/i 2 ) f 1 dx' f i 1 

kg + s 

[i?(k 1; k 2 , Q)$(a/, k', Q) - V(k u k 2 , Q)$(z', k, Q)] (5) 



^k, Q) ^ (, 2 ,g 2 ) + ^/;^/, 2 k^x 



where k = k' - k, k 1>2 = Q/2 ± k, k' lj2 = Q/2 ± k' and 

fl(k lf k 2 , Q) = + _ g 2 kg +^ (6) 

k'j + s k' 2 + s (k'i + s ) (k' 2 + s ) 



k' 2 + k 2 + 2s + k' 2 + k 2 + 2s ' 



v 0n, k 2 , Q) = r^-r^r- + t^tttt^T- (7) 



The following gluon propagator is used in ref. [|]: 

D < k > - < 8 > 

based on the original paper of Balitskii and Lipatov |3j , where the parameter s (which plays 
the role of a gluon mass) is generated using the Higgs mechanism. In [||, the parameter 
So is used to measure the magnitude of the infrared contributions to the cross sections. To 
scan these contributions, we will consider the following procedure: instead of the full BFKL 
Pomeron, we take the two gluon exchange and replace the perturbative gluon propagator 
by a non-perturbative one which includes the infrared contributions. 

The choice of the non-perturbative gluon propagators can be justified by the fact that 
they include infrared effects and then, they can be used in a more natural way to describe 
processes in which the infrared contributions play an important role. 



Before the use of the non-perturbative propagators, we will resume the main features of 
the methods to calculate the gluon propagator and discuss its distinct forms. 

The research of the infrared safe form of the gluon propagator has a long history ]13[ . The 
main trouble with the usual perturbative propagator is its divergence in the infrared region, 
more specifically, its pole at k 2 = 0. Them, obtain an infrared safe functional form for the 
gluon propagator several methods, both numerical and analytical are used: Dyson-Schwinger 
equations, lattice field theory and renormalization group methods. These frameworks give 
a number of solutions with different behaviors in the infrared and ultraviolet regions, which 
depend on the gauge considered and the approximations selected, and that we will discuss 
separately in the following. 

The Dyson-Schwinger equations (DSE) || are an infinite tower of coupled integral equa- 
tions among the Green functions of a field theory. For example, in QCD, the DSE for the 
quark self-energy involves the DSE for the quark and gluon propagators and the DSE for 
the quark vertex, and so on. In order to solve this system of equations, we must truncate it 
at some number of external legs and make an ansatz for the omitted part of the equations 
(normally using a Slavnov- Taylor identity). For example, in the case of DSE for the gluon 
propagator, we can restrict the number of external lines in the Green functions to three 
lines (no four gluon vertices), refuse the ghost propagator contributions and use the usual 
quark-gluon and three gluon vertices, resulting in a system with simplified equations for the 
gluon and the quark propagators. The solutions found in the literature (for the gluonic case) 
are very sensitive to the choice of the ansatz, to the choice of the gauge (Landau, Feynman, 
axial) and to the method used to solve the reduced set of equations, giving different behav- 
iors in the infrared region: infinite, finite (zero and non-zero). Some examples of the gluon 
propagator obtained with this technique are summarized below. 

• Cornwall |n| used a gauge independent resummation of Feynman diagrams (known as 
pinch technique) to obtain a special set of DSE's, whose solution presents an effective 
gluon mass. In [TJJ its value was determined as m^ ff = 500±200 MeV. The main feature 
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of this solution is a dynamically generated mass term in the propagator. Another 
remarkable feature of this solution is the correct ultraviolet behavior, according to the 
renormalization group. 



Habel et al. |15| used an alternative approach to calculate the Green functions for 
QCD. In the perturbation theory, the DSE's are solved by iteration, generating a power 



series for the propagators and vertices. In |15| , the principle of perturbation theory 
is maintained but some ansatz in the functional form of non-perturbative expression 
of the Green functions are made, which results in a simplified set of DSE's. This set 
yields the following solution for the gluon propagator: 

where b is a parameter to be determined. Note that the propagator vanishes at k 2 = 
and is finite when k 2 — > oo. Similar behavior was found by Gribov and Zwanziger [16 
using another arguments. 



Gorbar and Natale |L7j calculated the QCD vacuum energy as a function of the dynam- 
ical quark and gluon propagators using the effective potential for composite operators 
and using the operator product expansion (OPE) to relate the gluon and quark prop- 
agator with respectives condensates. The OPE gives the high energy behavior of the 
gluon polarization tensor II, which is related with the gluon propagator through (in 
the Landau gauge) 

v ^=-v^m(^-iP)- < 10 > 



In [|17| the authors use the following ansatz for the gluon tensor polarization (the 



momentum is in the Euclidean space), 



A 



= 0(X» 2 9 -K 2 ) + ^- 2 6(K 2 - W 2 ), (11) 



where 



and N is the number of colors (N = 3, in this case), ((cXs/^G^G^) ~ 0.01 GeV 4 is 
the gluon condensate |]TB||, resulting in fi g = 0.61149 GeV. The variational parameter 



X is determined numerically in [IT]] as x = 0.966797. 

We are aware of the use of the Lattice Field Theory (LFT) (for instance, see [0) as a non- 
perturbative method based on first principles, to make predictions about the infrared sector 
of QCD. However, LFT is not completely free of approximations. A lattice is characterized 
by two parameters: its size (number of points) and the separation between the points, and, 
due to the computational limitations, it is impossible to simulate an infinite lattice with an 
infinitesimal separation. Then a finite (and small) lattice with a given separation is used, 
introducing the called finite size and separation effects. Additional problems are the fermion 
duplication, solved by the quenched approximation and the phenomena of Gribov copies. 
Similarly to the DSE approach, the final result is also gauge dependent. In a recent paper, 



Alexandrou et al. [If] studied the gluon propagator with the lattice method in the Laplacian 
gauge - which is free of Gribov copies - and found that the best fit to the numerical results 
is obtained with the Cornwall's propagator. 

Having presented the features of the non-perturbative gluon propagators, we will employ 
them in a modified model to describe the process 77 — > J/^/J/^/ in the infrared region, using 



the propagators obtained using the DSE approach, specially the Gorbar-Natale one |17| and 



the Habel one [|Lj| 



As seen above, there are several different forms for the gluon propagator in the infrared 
region, nevertheless, these propagators can be tested in phenomenological applications. 
The most current model to describe the diffractive processes through non-perturbative 



gluons is the Landshoff-Nachtmann one |19[. In this model, the Pomeron consists of two 
non-perturbative gluons, whose usual properties are modified by QCD vacuum effects. This 
description of the Pomeron is very simple, since we don't include the multi-gluon exchange 
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effects, described by the BFKL equation ||. A consequence of this simplicity is the inde- 
pendence of the total cross-sections on the energy. Another problem is the fact that the 
model is based on an Abelian gauge theory, contrary to the non-abelian character of the 
QCD. Notwithstanding this, the LN model for the Pomeron, in connection with the Corn- 
wall's propagator, was used to describe the elastic pp scattering || and elastic production 
of vector mesons || with good agreement with experimental results. 

Here we use the model of two non-perturbative gluons for the Pomeron for the process 
77 — > J/^J/ty, and we consider the following approximations: first, we employ the Eq.(||) 
to identify the gluon propagator in the amplitude (Eq. @), and then we approximate the 
Pomeron as two non-perturbative gluons. With these assumptions, the amplitude is given 
by: 

(2)^3m#(M)= /— f$o(^,g 2 )l 2 ^kk-Q/2) 2 lp[(k+Q/2) 2 l (13) 

where T>(k 2 ) is the non-perturbative propagator. 

First, we substitute in the Eq.(|T3"D the Habel propagator (Eq. and vary the parameter 
b. This parameter is not determined in the works that proposed this propagator and we 
choose freely the values that give a result near to the one obtained in ||. Later, the results 
are compared with the propagator used in ||] (Eq. §), including the pure perturbative case 
(corresponding to sq = 0) as well as the pure massive case with a gluon mass similar to the 
one found by Cornwall [[L4[], So = 0.6 GeV 2 , since, in principle, we have freedom to set any 
value for it. As it can be seen in the results displayed in the fig.(U), the differential cross- 
section is sensitive to the choice of the propagator, but the global behavior is quite similar. 
The decrease of the cross section with both modified propagators is one order of magnitude or 
less, depending of the parameter involved and its value. When we compare the result of the 
perturbative case with the result of Eq.(|8|) with parameter so = 0.6 GeV 2 in \t\ = 5.0 GeV 2 , 
the cross-section decreases five times while for the Habel propagator with b = 0.6 GeV 2 at the 
same \t\ value, the cross-section is two times smaller. The overall behavior is a significative 
decrease in the cross-section with both propagators. A remarkable feature is the similarity of 



the results from different propag ators when s = 0.1 GeV 2 , b = 0.6 GeV 2 and s = 0.6 GeV 2 , 
b = 1.0 GeV 2 . The result is also compatible with the fact that the Habel propagator tends 
to the perturbative one when the parameter b goes to zero. 

In the next step we use the same strategy as above in the case of the Gorbar-Natale 
propagator, Eq. (|T0"D , but we can not vary freely the parameters since that \x s and x depend 
on the value of the gluon condensate, experimentally obtained. For this reason, we use the 
values for fi g and \ calculated in [D]]. The result is shown in the fig.(f|) in comparison with 
the perturbative propagator, and is quite similar to the Habel propagator: with the choice of 
the parameters above the non-perturbative propagator gives a result near the perturbative 
one with a massive term with sq = 0.1 GeV 2 . When this non-perturbative propagator is 
employed, the decrease of the cross-section in \t\ = 0.5 GeV 2 is around 40 %. 



In the previous works in which the non-perturbative gluon propagators are used the 
running coupling constant is frozen at some value when Q 2 — > 0. In the above results, we use 
a running coupling constant with a k 2 dependence when Q 2 goes to zero. In order to test the 
sensitivity of the cross-section with respect the choice of the coupling constant we computed 
the same observable, but now with a coupling constant with the above requirement, following 



the result from Cornwall 114 



*s(Q 2 ) = (14) 

251n(^^) 

for the case of the Gorbar-Natale propagator. In the fig.© we shown the results, where 
the change of the coupling constant increases considerably the differential cross section 
(approximately two orders of magnitude) and modifies crucially the slope in the low values 
of t. 

The results obtained show that the behavior of the cross-section could be quite different 
(we obtain a decrease of one order of magnitude) once we change the gluon propagator. 
However, without experimental results expected in TESLA, we cannot restrict reasonably the 
choice of the propagator and the parameters for this process. The comparison with another 
diffractive process calculation obtained through the same methods could be employed, for 
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example the elastic pp(pp) scattering or elastic photo-production of vector mesons, which 
we will soon consider. 

Someone can argue against the use of the non-perturbative propagator in a perturba- 
tive calculation, but we only consider the first term in the perturbative expansion and the 
calculation of the form factor for the transition photon-vector meson is calculated in a non- 
perturbative manner, therefore our calculation is valid in the kinematical region of interest. 
More clearly, we are only interested in an estimation of the effects in the infrared region. We 
can see that with a particular choice of the parameters, the result with a bare gluon mass @] 
and the another one, with a non-perturbative propagator, is very similar (see figs.@ and 
(f|) ) . It is well known that a massive gauge field theory like massive Yang-Mills presents sev- 
eral problems (see 0] for a review and use of this theory to describe diffractive phenomena) 
as well as the non-perturbative propagators but in the infrared region the last one seems to 
be a better choice than the former to describe diffractive processes. 

From our results, we can see that these effects in the region \t\ < 2 GeV 2 diminish the 
differential cross-section for the distinct propagators once compared with the pure perturba- 
tive result. The goal of this work is to analyze the effects in the cross-section of the process 
77 — > J/^/J/^ in the infrared region induced with the replacement of the usual gluon prop- 
agator by a distinct one, obtained from non-perturbative methods, which is more consistent 
in the kinematical region of interest. We show that the differential cross-section is reduced 
in comparison with the pure perturbative result in the very low momentum transfer regime 
with a running coupling constant. However, we should only be able to select which is the 
propagator that describes better the process once experimental results are provided. 
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FIGURES 
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FIG. 3. The differential cross-section of the process 77 — > J/^f J/\l/ using the gluon propagators 
of Lipatov et al. (Eq.(|8|), parameter sq) [f§ and Habel et a/.(Eq.@, parameter b). 
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FIG. 4. The differential cross-section of the process 77 — > J/^f J/ty using the gluon propagators 
of Lipatov et al. (Eq.(||), parameter sq) @ and Gorbar and Natale (Eq.(|i~0"l)) with the parameters 
obtained in fp~7p| . 
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FIG. 5. Comparison between the results for differential cross-section obtained with a running 



coupling constant and frozen coupling with the Gorbar-Natale propagator [17| 
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